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Starting from a microscopic lattice model, we investigate clustering, micellization and micelle 
ordering in semi-dilute solutions of AB diblock copolymers in a selective solvent. To bridge the 
gap in length scales, from monomers to ordered micellar structures, we implement a two-step coarse 
graining strategy, whereby the AB copolymers are mapped onto "ultrasoft" dumbells with monomer- 
averaged effective interactions between the centres of mass of the blocks. Monte Carlo simulations 
of this coarse-grained model yield clear-cut evidence for self-assembly into micelles with a mean 
aggregation number n ~ 100 beyond a critical concentration. At a slightly higher concentration 
the micelles spontaneously undergo a disorder-order transition to a cubic phase. We determine the 
effective potential between these micelles from first principles. 



Block copolymers in solution show a remarkable ten- 
dency to self-assemble into a bewildering range of dis- 
ordered (liquid-like) and ordered structures, depending 
on the macromolecular composition, the relative sizes of 
the blocks, solvent selectivity, polymer concentration and 
temperature 0, Q . Complex phase diagrams of styrene- 
isoprene diblock copolymers in various solvents have been 
mapped out experimentally in great detail by varying the 
above physical parameters, using a combination of light 
scattering, small-angle neutron and X-ray scattering and 
rheological studies [3|. A frequent, partial scenario sees 
the copolymers aggregate into polydisperse spherical mi- 
celles at fairly low polymer concentration; upon lowering 
the temperature or increasing the concentration, the mi- 
celles undergo a disorder-order transition onto a cubic 
lattice. A theoretical understanding of copolymer phase 
behaviour generally relies on phenomenological consider- 
ations , or on self-consistent field theory similar to that 
applied earlier to copolymer melts 6] . A more molecular 
approach is obviously very difficult, in view of the wide 
range of length scales involved, from the monomer level 
to the mesoscopic scales characterising ordered micellar 
structures. In particular, monomer-level Monte Carlo 
(MC) simulations have so far been restricted to short 
copolymer chains 0. 

In this letter we propose and implement a two-stage 
coarse graining strategy to investigate the micellization 
and disorder-order transition of a simple microscopic 
model of diblock copolymers. The AB copolymers are 
made up of two blocks of equal numbers Ma = Mb = M 
of monomers on a cubic lattice. N such chains occupy a 
fraction of the lattice sites. The remaining sites are oc- 
cupied by an implicit solvent chosen to be selective (good 
solvent) for the B blocks, while it is a ^-solvent for the 
A blocks. This situation is crudely represented by mod- 
elling the A blocks as ideal chains (I) and the B blocks as 
self and mutually avoiding walks (S). Moreover A and B 
blocks of the same or different copolymers are assumed 



to be mutually avoiding (S). The resulting model will be 
referred to as ISS. The model is clearly athermal, so that 
temperature will be irrelevant, and the polymer density 
p = N / L 3 is the only thermodynamic variable (L 3 being 
the total number of sites of the cubic lattice, the spacing 
a of which will be taken as the unit of length) . With typ- 
ical polymer sizes M ~ 10 3 , and expected micelle aggre- 
gation numbers n ~ 10 2 , very large system sizes would 
be needed to generate the tens of micelles required to 
study their ordered and disordered structures. 

The first stage of our coarse-graining procedure, de- 
signed to overcome the above bottleneck, is inspired 
by earlier work on homopolymers [8|, |9j. AB copoly- 
mers are mapped onto "ultrasoft" asymmetric dumbells 
with effective interactions v a p{R) between the centres of 
mass (CM) of A and B blocks on different copolymers, 
and an intramolecular "tethering" potential 4>ab {R) be- 
tween the CM's of the two blocks of the same copoly- 
mer Since the original model is athermal, these 
effective interactions are purely entropic and are deter- 
mined by taking statistical averages over copolymer con- 
formations for fixed A-A, A-B and B-B CM distances 
R. In practice this is achieved by inverting the CM- 
CM intermolecular and intramolecular pair distribution 
functions gAA{R), 9ab{R), 9bb{R) and sab(R)- The in- 
version was carried out in the low density limit, using 
pair distribution functions generated in monomer-level 
MC simulations of two diblock copolymers with 2M = 
500 monomers each, and an exact relation between the 
Qnfi(R), sab(R) and the corresponding v a p(R),<j>AB(R) 
[lOj. The resulting pair potentials are shown in figure fj] 
where they are compared to the corresponding pair po- 
tentials for an equimolar binary mixture of untethered A 
and B homopolymers; in the latter case the inversion pro- 
cedure is trivial in the zero density limit where the effec- 
tive potentials v a p(R) and 4>ab{R) reduce to the poten- 
tials of mean force, e.g. v a p(R) = —kBT\ng a p{R). The 
two sets of effective pair potentials vab(R) and vbb(R) 
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FIG. 1: (Colour on line) Effective pair potentials 
vaa{R),vab{R),vbb{R) and 4>ab{R) for tethered and un- 
tethered block copolymers in the zero density limit. 



are seen to be surprisingly close. As in the case of ho- 
mopolymers 0, |tj they are roughly gaussian in shape 
with an amplitude of a few ksT, and a range of the order 
of the natural length scale, namely the radius of gyration 
R g of the copolymer. In lattice units, the MC simulations 
yield the following radii for an isolated copolymer with 



2M = 500: R 



g,A 



6.964, R„ B = 11.274, R„ 



13.889, 



while the root mean square distance between the CM's of 
the two blocks is Rab = 20.661. The results for vab(R) 
and Vbb(R) in figure Q] suggest that tethering of the two 
blocks does not significantly modify their effective inter- 
actions. This is no longer true for vaa(R), since this po- 
tential is strictly zero in the untethered (binary mixture) 
case, while it is about 10% of vbb(R) for the copolymer; 
this residual effective interaction between the ideal blocks 
is induced by the tethering to non-ideal B-blocks. 

The key assumption made in this work is that these 
zero density pair potentials are transferrable to finite 
copolymer densities, i.e. we neglect the density depen- 
dence of the effective potentials. Earlier experience with 
homopolymers shows that this is not an unreasonable 
approximation well into the semi-dilute regime p/p* 5 
where p* = (47i\R^/3) -1 is the overlap density 0. 

We have carried out extensive MC simulations of sam- 
ples of 5 x 10 3 and 10 4 AB copolymers using the effective 
pair potentials of figure ^ in the range 3 < p/p* < 10. 
The simulations using the coarse-grained representation 
are at least two order of magnitude faster than simu- 
lations of the original all-monomer lattice model. This 
allowed us to make detailed investigations, using both 
sets of potentials in figure in order to ascertain the 
sensitivity of the micellization behaviour to details of the 
effective pair potentials. Both sets of potentials yield re- 
sults in semi-quantitative agreement, and only those ob- 
tained with the untethered potentials will be shown here. 
Inspection of snapshots of configurations shows that mi- 



cellization sets in between p/p* = 4 and p/p* — 5, with 
the ideal A-blocks forming the dense core, as one might 
expect because of the low or vanishing energy penalty, 
while the B-blocks form the corona, where their energet- 
ically unfavourable overlap is minimized. An unambigu- 
ous definition of a micelle remains somewhat arbitrary 
because of the size polydispersity, but we used the op- 
erational convention that all A-blocks whose centres are 
sit 
to 
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FIG. 2: (Colour on line) Panel a): probability of a cluster of 
n monomers. Panel b): probability distribution of the cluster 
radius R c - 

aggregation numbers n are shown in panel a of figure|21for 
p/p* = 3, 5 and 6, based on a cut-off distance r c = 0.8R g . 
P(n) turns out to be very sensitive to the precise value 
of r c for p/p* < 4 but nearly independent of the cut-off 
at higher densities, suggesting that well defined micelles 
appear only beyond an estimated cmc of p/p* ~ 5. At 
that concentration P(n) exhibits a broad peak around 
n ~ 70, which narrows and shifts to a mean aggregation 
number n ~ 100 for p/p* > 6. Panel b of figure [3 shows 
the corresponding distribution P(R C ) of the cluster radii 
(R c ) defined as 
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R: 



Rcm 



(1) 



where R A is the position of the i-th A-block belonging to 
an n-cluster and Rcm the position of the instantaneous 
CM of that cluster. At p/p* — 3, P(R C ) is seen to peak 
at a low value of R c corresponding to small clusters; the 
distribution becomes bimodal at p/p* — 5, signalling the 
proximity of the cmc, while at p/p* = 6, P{R C ) shows a 
sharp peak around R c ~ R gi corresponding to micelles 
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of well defined core size. The peak sharpens and moves 
to lower values of R c at still larger values of p/p*, in- 
dicative of a contraction of the core with copolymer con- 
centration. The internal structure of the micelles is best 
( if 
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FIG. 3: p/p* — 7. Density profiles of A-blocks (open sym- 
bols) and B-blocks (closed symbols) with respect to the CM 
of each cluster averaged over all clusters with n > 40. Note 
the different scales for pA (left axis) and pb (right axis). 

over all micelles formed during the MC simulations, i.e. 
over size polydispersity. As expected, pa(R) peaks at 
the origin and drops rapidly to zero for R > R g , while 
the profile pb{R) practically vanishes for r < R g and 
peaks at R ~ 2R g , thus defining the corona radius. The 
sharp distinction between micelle core and corona evi- 
dent from profiles at p/p* = 7 gradually blurs as the 
copolymer concentration decreases thus signalling that 
the micelles dissolve into a quasi-homogeneous solution 
of copolymers. 

A periodic sample of 5000 AB copolymers generated 
about 50 micelles, and the second step of our coarse- 
graining procedure is to consider the structure of the 
system of micelles, regarded as well defined spherical 
colloids. Figure 01 shows the pair distribution function 
gmm(R) of the CM's of the micelles, extracted from the 
MC simulations for copolymers concentrations p/p* = 
5, 5.5 and 6. All g mm {R) go rapidly to zero for R < 3R g , 
indicative of a strong effective short-range repulsion be- 
tween micelles. The structure at p/p* = 5 is typical of 
that of a dense colloidal fluid, but it is considerably en- 
hanced at p/p* = 6, with the amplitude of the main peak 
jumping from 1.7 (at p/p* = 5) to over 3 while the split 
second peak in g m m{R) is characteristic of a cubic struc- 
ture. The sudden change in pair structure is suggestive 
of a disorder /order transition, and snapshots of micelle 
configurations, shown in figure 01 confirm this interpre- 
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FIG. 4: (Color on line) Pair distribution function g m m{R) 
of the CM's of the micelles at p/p* = 5 (closed circles), 5.5 
(open squares) and 6 (closed squares). The HNC potential 
Vm—m(R) at p/p* = 5 is shown by a continuous line (in red). 
Also the corresponding pair correlation is shown as a thick 
curve in very good agreement with g m m(R) at p/p* — 5. 
Note that the potential scale in on the right axis. Finally, the 
snapshots in the upper part of the figure show the x-y plane 
projection of the position of the centre of mass of the micelles 
in a typical configuration at p/p* = 5 (left-disordered) and 6 
(right-ordered). 

tation. Note that at p/p* = 6, an initially disordered 
micelle configuration evolves spontaneously to the defec- 
tive ordered structure shown in figure Rafter few tens of 
thousands MC cycles. 

From the MC generated pair distribution function 
9mm{R) at p/p* = 5 one can extract an effective micelle- 
micelle pair potential using the HNC integral equation of 
liquid state theory [n], according to 

v mm (R) = k B T {g mm (R) - 1 - c mm (R) - \n[g mm {R)]} 

(2) 

where the direct correlation function c mm {R) is related 
to g m m(R) by the Ornstein-Zernike relation The re- 
sulting effective pair potential shown in figure 0] exhibits 
a relatively soft repulsion, which culminates at a finite 
value at r = 0, followed by a shallow attractive well. A 
precise value for the potential at the origin is difficult to 
obtain because of the noise in g mm at short distance. As- 
suming a gaussian-like core, a fit to the data in the range 
0.5 < r < 3.5 provides v mm (0) ~ \2ksT. The accuracy 
of the inversion procedure is vindicated by a direct MC 
simulation of micelles (represented as point particles) in- 
teracting via this v mm (R): the resulting pair distribution 
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functions, also shown in figure^] is indistinguishable from 
the g mm (R) generated in the original simulation of the 
underlying copolymer system, signalling the success of 
the second step in our hierarchical coarse-graining pro- 
cedure. Potentials like v mm (R) shown in figure 4, which 
remain finite as R — > 0, can give rise to a freezing tran- 
sition at sufficiently high densities, as shown, e.g., in the 
case of the gaussian core potential 0, 0|. The over- 
all scenario described above is also observed when the 
low density effective potentials v a p(R) for the tethered 
AB copolymer are used (rather than the untethered bi- 
nary mixture potential, cfg. figure [IJ, except that the 
cmc and disorder-order transition are shifted to higher 
densities (p/p* — 6 and 7 respectively). Moreover the 
A-block density profiles (similar to those shown in figure 
EJ, now develop a minimum at R = at higher densi- 
ties, so that the micelles tend to become "hollow" and 
increase in size. This is easily understood from the weak 
short-range repulsion in Vaa (R) shown in figure ^ 

In summary, we have implemented a two-step coarse- 
graining procedure for the description of self-assembly 
and micelle ordering of a simple model for diblock copoly- 
mers in a selective solvent. In the first step, the diblock 
chains are mapped onto ultrasoft dumbells, and the inter 
and intramolecular potentials are determined by invert- 
ing low-density structural MC data for the initial lattice 
model. MC simulations of large samples of the coarse- 
grained dumbell model above the overlap density show 
cluster formation, micellization immediately followed by 
ordering of the micelles onto a cubic lattice, in qualita- 
tive agreement with experimental observation 
In a second step we have determined the effective in- 
teraction between the moderately polydisperse spheri- 
cal micelles, which may be put to good use in a more 
quantitative investigation of micelle ordering. Our coarse 
graining methodology raises the general question of the 
transferability of the low-density effective potentials to 
semi-dilute copolymer solutions, which we are presently 
investigating. The present athermal model of a diblock 
copolymer is clearly an over-simplification, but tempera- 
ture effects may be easily included by introducing interac- 
tions between adjacent monomers on the lattice IJ. The 
related symmetric ISI model investigated in 01 > where A 
and B are ideal, but mutually avoiding, exhibits a lamel- 
lar phase rather than the present micellar phase, illus- 
trating the strong influence of solvent selectivity. So far 
we have considered block of equal lengths but the gener- 
alisation of our methodology to different relative lengths 
/ of the blocks is straightforward, and is being pursued. 
Finally we believe that coarse-graining provides new in- 



sight into clustering and micellization, as a consequence 
of the competition between entropy and low energy over- 
lap penalties [T H Hr l Il7 t . 
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